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We formulate equations of time-dependent density functional theory (TDDFT) in the co-moving 
Lagrangian reference frame. The main advantage of the Lagrangian description of many-body dy- 
namics is that in the co-moving frame the current density vanishes, while the density of particles 
becomes independent of time. Therefore a co-moving observer will see the picture which is very 
similar to that seen in the equilibrium system from the laboratory frame. It is shown that the most 
natural set of basic variables in TDDFT includes the Lagrangian coordinate, £, a symmetric defor- 
mation tensor g^, and a skew-symmetric vorticity tensor, Fp V . These three quantities, respectively, 
describe the translation, deformation, and the rotation of an infinitesimal fluid element. Reformula- 
tion of TDDFT in terms of new basic variables resolves the problem of nonlocality and thus allows 
to regularly derive a local nonadiabatic approximation for exchange correlation (xc) potential. Sta- 
tionarity of the density in the co-moving frame makes the derivation to a large extent similar to the 
derivation of the standard static local density approximation. We present a few explicit examples 
of nonlinear nonadiabatic xc functionals in a form convenient for practical applications. 
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I. INTRODUCTION 

Traditionally, the physical understanding of various 
many-body phenomena is based on Landau's intuitive 
concept of quasiparticles, which relates the behavior of 
a strongly interacting quantum system to the properties 
of a gas of noninteracting (or weakly interacting) elemen- 
tary excitations 1 . The common field-theoretical formula- 
tion of the many-body problem 2,3,4 allows to rigorously 
justify this very appealing point of view, provided the 
system is in a weakly exited many-body state. Unfor- 
tunately, practical applications of traditional many-body 
methods to real systems, even in the equilibrium or lin- 
ear response regime, is computationally very demanding. 
In strongly nonequilibrium systems the situation is much 
worse. The simple intuitive picture of elementary exci- 
tations breaks down, while the direct application of the 
standard many-body theory becomes increasingly diffi- 
cult even for model systems. 

An alternative approach to the quantum many-body 
problem is offered by the density functional theory 
(DFT) 5 ' 6 ' 7 . DFT represents a mathematically rigorous 
realization of another famous idea in theoretical physics 
which is a concept of collective variables theory 8 . In- 
deed, DFT opens a possibility to formulate the many- 
body problem in a form of a closed theory that contains 
only a restricted set of basic variables, such as density 
in the static DFT, 5 or density and current in the time- 
dependent DFT (TDDFT) 6,7,9 . In classical physics a the- 
ory of this type is known for more than 250 years. This is 
classical hydrodynamics. In fact, the Runge-Gross map- 
ping theorem in TDDFT 9 proves the existence of an ex- 
act quantum hydrodynamics. An analogy of TDDFT to 
hydrodynamics has been already noted in the original 
paper by Runge and Gross, Ref. 9 (see also our recent 



paper, Ref. 10). In this respect the static DFT should 
be viewed as an exact quantum hydrostatics. It is in- 
deed known that the condition of energy minimum is 
equivalent to the condition for a local compensation of 
the external and internal stress forces, exerted on every 
infinitesimal element of an equilibrium system 11 . Thus 
DFT not only introduces an alternative formalism in the 
quantum many-body problem, but it also naturally sug- 
gests an alternative way of thinking, which refers to the 
physical intuition, developed over hundreds year of ex- 
perience in the classical continuum mechanics. Interest- 
ingly, equations of TDDFT in the hydrodynamic formu- 
lation can be also considered as a force balance condi- 
tion, but in a local noninertial reference frame moving 
with the flow. In the time-dependent case there is a lo- 
cal compensation of the external, inertial, and internal 
stress forces. This demonstrates a close similarity of the 
static DFT (which is currently a well developed theory) 
and TDDFT (which is still under construction) in the co- 
moving frame. The above similarity was the main mo- 
tivation to reconsider the formulation of TDDFT from 
the point of view of a local observer in the co-moving 
Lagrangian reference frame. In this paper we present 
the results of such a reconsideration, based on our recent 
formulation of the many-body theory in the Lagrangian 
frame 12 (in what follows the papers of Ref. 12 is referred 
to as I). 

Practical applications of any DFT rely on the Kohn- 
Sham (KS) construction 5,13 , which maps the calculation 
of basic observables in the interacting system to the solu- 
tion of an auxiliary, noninteracting KS problem. Nonin- 
teracting KS particles move in a self-consistent exchange 
correlation (xc) potential that is adjusted to reproduce 
the correct values of basic variables, i. e. density and 
current in TDDFT. From the hydrodynamical point of 
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view the KS construction allows one to compute exactly 
the kinetic part of the internal stress force, while treat- 
ing the xc contribution to the stress in an approximate 
fashion. Thus the central problem of any practical DFT 
reduces to the construction of adequate approximations 
for xc potentials. In the static DFT a good starting point 
is provided by the local density approximation (LDA). 
On the one hand, the static LDA, by itself, gives quite 
reasonable results, and, on the other hand, it allows for 
further modifications and refinements. The construction 
of a similar basic local approximation in TDDFT turned 
out to be problematic. The reason for these problems 
is the inherent nonlocality of the nonadiabatic nonequi- 
librium DFT. General arguments 14 ' 15 , based on the har- 
monic potential theorem 16 , require that any consistent 
nonadiabatic xc potential must be a strongly nonlocal 
functional of the density. Otherwise the harmonic po- 
tential theorem is violated. 

In 1996 Vignale and Kohn 17 (VK) realized that at least 
in the linear response regime the problem of nonlocality 
can be resolved by changing the basic variable and by 
using the xc vector potential, A xc , instead of the scalar 
one. Namely, VK showed that if one considers the cur- 
rent j (instead of the density n) as a basic variable, a 
consistent linear local approximation for A xc can be reg- 
ularly derived 17 . Shortly after that, Vignale, Ullrich and 
Conti 18 (VUC) found that the further switching of vari- 
ables from the current j to the velocity v = j/n, al- 
lows to represent the complicated VK expression in a 
physically transparent viscoelastic form. VUC also phe- 
nomenologically constructed a nonlinear nonadiabatic xc 
functional. In this construction they formally adopted 
the linear VK form, but with coefficients, taken at a "re- 
tarded position" 19 . A similar, linear in velocity gradi- 
ents, phenomenological construction based on Keldysh- 
contour action functional, has been recently proposed by 
Kurzweil and Baer 20 . An attempt to regularly derive a 
nonlinear nonadiabatic local approximation in TDDFT 
was made in our recent work 10 . Noting that the applica- 
bility conditions for the linear VK approximation exactly 
coincide with those for the "collisionless hydrodynamics" 
of Rcfs. 21,22,23, we extended the hydrodynamics formal- 
ism to TDDFT. Using Landau Fermi-liquid theory we 
derived a nonadiabatic nonlinear approximation for the 
xc stress tensor that defines the xc potential 10 . It has 
been shown that the stress tensor is a local functional 
of new basic variables: the Lagrangian coordinate and 
a second-rank metric-like tensor. An appearance of the 
Lagrangian coordinate as a basic variable is quite remark- 
able if we recall the abovementioned (see also I) static-like 
force balance in the co-moving Lagrangian frame. 

In this paper we reformulate TDDFT using the exact 
equations of quantum many-body dynamics in the La- 
grangian frame 12 . We show that possibly the most natu- 
ral complete set of basic variables in TDDFT consists of 
the Lagrangian coordinate £, the symmetric Green's de- 
formation tensor g^ Vl and the skew-symmetric vorticity 

tensor F^ u . These three quantities, one vector, one sym- 



metric and one skew-symmetric tensors contain twelve 
numbers that are required for the complete local charac- 
terization of a deformed state of any continuum media 24 . 
Namely, £, g^ v and F^ respectively describe the trans- 
lation, deformation and the rotation of an infinitesimal 
fluid element. On the other hand, tensors g^ v and F^ v de- 
scribe generalized inertia forces in the Lagrangian frame 
(see I). All three quantities are functionals of velocity 
in accordance with the Runge- Gross mapping theorem 9 . 
The new formulation of TDDFT relates the local stress 
in the system to the dynamic deformations, which is very 
natural physically. The main practical advantage is that 
the reformulation of TDDFT in terms of new basic vari- 
ables resolves the problem of nonlocality on the most 
fundamental level. We show that the dynamic xc stress- 
deformation relation can be made local. This allows us 
to derive a local nonadiabatic approximation in a reg- 
ular way that is similar to the derivation of the static 
LDA. The whole history of nonadiabatic approximations 
in TDDFT can be viewed as a staircase of successive 
transformations of basic variables, from the density, via 
the current and the velocity to the general geometric 
characteristics of deformed media. The first step was 
made by VK 17 in 1996. Hopefully the present formu- 
lation of TDDFT corresponds to the last step on this 
staircase. 

The structure of the paper is the following. In Sec. II 
we consider the hydrodynamic formulation of TDDFT. 
Using this formulation we introduce the KS system and 
define xc potentials in terms of stress tensors. In Sec. Ill 
we discuss local approximations in DFT. We derive the 
common static LDA and explain physical reasons of the 
nonlocality in TDDFT. Sections IV and V are the cen- 
tral parts of the present paper. In Sec. IV we introduce 
new basic variables and develop a complete geometric 
formulation of TDDFT. We also formulate a universal 
many-body problem which determines the xc stress ten- 
sor and the xc potential as functionals of basic variables. 
In Sec. V we discuss an approximate solution of the uni- 
versal many-body problem in the lowest order of the gra- 
dient expansion. Explicit examples of nonadiabatic local 
xc functionals are presented in Sees. VB-VD. In the con- 
cluding part of Sec. VD we formulate the full set of time- 
dependent KS equations in the nonlinear time-dependent 
local deformation approximation. In equilibrium case 
these equations reduce to the common equations of DFT 
with the static LDA xc potential, while in the linear re- 
sponse regime we recover the results of VK approxima- 
tion. In Sec. VI we summarize our results. Three ap- 
pendixes contain technical details of calculations. 



II. DEFINITION OF EXCHANGE 
CORRELATION POTENTIALS VIA STRESS 
DENSITY 

In this section we discuss a hydrodynamic formulation 
of DFT and introduce a definition of xc potentials in 
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terms of local stress forces. Let us consider a system of N 
interacting fcrmions in the presence of a time-dependent 
external potential £/ ex t (x, t) . The system is described by 
the following Hamiltonian 



H = T + W + U, 



(1) 



where T, W, and U correspond to the operators of the 
kinetic energy, interparticle interaction energy, and the 
energy of interaction with the external field respectively 

T = -yW(x)|-V>(x), (2) 

W=^ J dxdx'w(\x-x'\)^(x)^(x')ip(x')il>(x), (3) 

U= f dxU ext {x,t)^(x)^{x). (4) 

Here w(|x|) is the potential of a pairwise interparticle in- 
teraction. Field operators ip^ and ip satisfy the fermionic 
commutation relation 



{^( X )^(x')} = <5(x-x'). 



(5) 



For definiteness we consider a Fermi system, although 
most of the results of this paper are independent of statis- 
tics. 

The hydrodynamic formulation of DFT follows the 
Heisenberg equations of motion for the density and for 
the current operators 



^-n(x,t)-i([H,n(x)})=0, 
|j(x,i)-i([tfj(x)])=0, 



(6) 
(7) 



where n(x) = (n(x)) and j(x) = (j(x)) arc the density 
of particles and the current density. The corresponding 
operators are defined by the standard expressions 

n{x,t)=^{x,t)<P(x,t), (8) 

Angle brackets in the above formulas denote averaging 
with the exact density matrix jo, i.e. (A) — TrpA. Equa- 
tions (6), (7) can be represented in a form of the hydro- 
dynamics balance equations (details of the derivation can 
be found, for example, in I) 



d 

Dtn + n d^ V » =0 > 



(10) 



mnDfVfj + — P, v + n—U = 0, (11) 



Equation (10) is the common continuity equation, while 
Eq. (11) corresponds to the local momentum conserva- 
tion law. In these equations v = j/n is the velocity of 



U = U cxt + Uh_ is the sum of the external and the Hartree 
potentials, 

U H (x,t) = J w(\x-x'\)n{x',t)dx'. (12) 

The exact stress tensor in Eq. (11), 

(x,t)=T„ v (x,t) + W llv (x,t), (13) 

contains the kinetic, T^ v , and the interaction, W^, con- 
tributions. Divergences of the tensors T^ u and in 
Eq. (11) come from the commutators of the current op- 
erator with operators T, Eq. (2), and W, Eq. (3), respec- 
tively. In I we have derived the following explicit repre- 
sentations for the stress tensors (see also Refs. 2,25,26,27) 

^(x) = ±-({q^q v ^ + (q^)%yj - ^V 2 n) (14) 

1 f ,xV dw(\x'\) 
W„„(x) = --Jdx |x/| 9|x1 

x f G 2 (x + Ax',x- (1- \)x')d\ (15) 
Jo 

where q = — iV — rav is the operator of "relative" mo- 
mentum, and G2(x, x') = (x)n(x')^(x)} — n(x)n(x') 
is the pair correlation function 28 . It is worth mentioning 
that the representability of the stress force in a form of 
a divergence of a tensor is a direct consequence of the 
Newton's third law (see Appendix A in I). 

Equations (10) and (11) represent the exact local 
conservation laws which must be satisfied for an arbi- 
trary evolution of the system. Let us apply them to 
TDDFT. The first, less restrictive part of TDDFT map- 
ping theorem 9,29 states the existence of a unique and in- 
vertible map: j — > U cxt or, equivalently, v — > C/ ext . This 
implies that the exact many-body density matrix p(t) for 
a given initial condition, p(0) = po, is a functional of the 
velocity v. Hence the stress tensor of Eqs. (13) is a func- 
tional of v and of the initial density matrix: P^ v [po , v] . 
Therefore Eqs. (10) and (11) constitute a formally closed 
set of the exact quantum hydrodynamics equations with 
the memory of initial many-body correlations. It is inter- 
esting to note that the common classical hydrodynamics 
can be viewed as a particular limiting form of TDDFT. 
In this limiting case the stress tensor functional is known 
explicitly - it takes the usual Navier-Stokes form 30 . 

In the equilibrium system Eq. (11) reduces to the static 
force balance equation 



_d 

dx 



d_ 



-P^+n^—U = Q 



(16) 



the flow, D t 



dt 



vV is the convectivc derivative and 



This equation shows that the force produced by the ex- 
ternal and the Hartree potentials is compensated by the 
force of internal stresses. The net force, exerted on every 
infinitesimal fluid clement is zero, which results in zero 
current density and a stationary particles' density distri- 
bution. According to the Hohcnberg-Kohn theorem 5 ' 31 
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any equilibrium observable, in particular the stress ten- 
sor, is a functional of the density n: P^ v = P^ u [n\. Hence 
Eq. (16) is, in fact, the equation of the exact quantum hy- 
drostatics that uniquely determines the density distribu- 
tion in the equilibrium system. In the semiclassical limit 
Eq. (16) reduces to the common hydrostatics equation 30 
(for a degenerate high density Fermi gas we recover the 
Thomas- Fermi theory). 

Practical application of DFT always rely on the KS 
construction that is a particular consequence of the map- 
ping theorems. The current and the density in the inter- 
acting system can be reproduced in a system of nonin- 
tcracting KS particles, moving in a properly chosen self- 
consistent potential Us = U + U xc 32 - The hydrodynamic 
formulation of TDDFT/DFT allows us to relate the xc 
potential U KC to the stress density. Hydrodynamics bal- 
ance equations for the KS system take the form 

^ + "^,=0, (17) 

mnDtV » + 3^ + H ^ Uxc + = ° (18) 



dx>* 



where the kinetic stress tensor of KS system, , is given 
by Eq. (14), but with the averaging over the state of 
nonintcracting particles. Comparing Eq. (10) and (11) 
with Eqs. (17) and (18) we find that the velocity v and 
the density n of the noninteracting and the interacting 
systems coincide if the xc potential t/ xc (x, t) satisfies the 
equation 



1 <9P XC 
dx^ n dx v 



(19) 



where P x £ is the xc stress tensor that equals to the dif- 
ference of the stress tensors in the interacting and non- 
interacting systems with the same velocity distribution 



pxc p rp^ 



(20) 



Equations (19) and (20) demonstrate the physical signif- 
icance of U xc . The xc potential should produce a force 
which compensates the difference of the internal stress 
forces in the real interacting system and in the auxil- 
iary noninteracting KS system. By continuity equation, 
the density n is a functional of the velocity. Therefore 
Eq. (19) defines U xc (up to inessential constant) as a func- 
tional of v. Equation (19) shows that n _1 <9„P x £ is a po- 
tential vector. This does not mean that vectors n^ 1 d„P flv 
and n~ 1 duT^ i/ separately have no rotational components. 
However, according to the balance equations of Eqs. (11) 
and (18) the rotational components of these vectors co- 
incide - both equal to the rotational part of the vector 
mD t v. 

It is also possible to construct the proper KS sys- 
tem using the xc vector potential A xc or, in the most 
general case, a combination of xc vector and scalar 
potentials 10,17 ' 18,33 . In this case the exact local conser- 
vation laws require that the total xc force, ^F xc , should 



compensate the difference of the stress forces in the in- 
teracting and the KS systems: 



■pxc 



dA x u c dU. 



1 <9P XC 
dx^ n dx v 



(21) 



This equation determines the xc potentials, A xc and C/ xc , 
up to a gauge transform. Equation (21) represents a very 
important exact property of xc potentials: they produce 
a force that must be a divergence of a second rank tensor. 
This requirement automatically implies the well known 
zero force and zero torque sum rules 6 



J n^ xc dx = 0, J n(x x 



)dx = 0. (22) 



We would like to outline that the exact local condition 
of Eq. (21) is much stronger then the common integral 
requirements of Eq. (22). Apparently the above defini- 
tion of xc potentials equally well apply both to TDDFT 
and to the static/equilibrium DFT. It should be men- 
tioned that in the equilibrium case the stress forces in 
the interacting and KS systems separately are potential 
vectors. 

Let us briefly discuss the role of xc vector potential in 
DFT. Apparently an appearance of A xc is unavoidable 
in the presence of an external magnetic field 33 . Indepen- 
dently of the character of external fields, the formula- 
tion in terms of A xc is convenient in the linear response 
regime 17,18 . Indeed, in the linearized theory we can per- 
form the Fourier transform in the time domain, which 
makes A xc completely local, provided the xc stress ten- 
sor is a local functional of some basic variables. In the 
nonlinear regime this advantage clearly disappears. For a 
nonlinear evolution the description of xc effects in terms 
of the scalar potential, defined by Eqs. (23), (24), looks 
at least as convenient as the formulation in terms of A xc . 
Below for definitencss we assume the noninteracting v- 
representability of the velocity, which allows us to con- 
struct the KS system using only the scalar xc potential. 
Reformulation of the theory for xc vector potential is 
straightforward (see Conclusion). 

For the practical applications it is possibly more con- 
venient to represent the force definition of U xc , Eq. (19), 
in a familiar form of the Poisson equation 



V 2 t/ xc (x,i) = 4 7 T ( O xc (x,f), 



where the quantity p xc (x, i), 

1 d ( 1 



Pxc = 



Air dx^ 



9 pxc 

ndx" » v 



(23) 



(24) 



can be interpreted as an xc "charge" density. In this 
context the xc stress force, n~ 1 d u P*£,, has a clear mean- 
ing of an xc "polarization" density. The additional dif- 
ferentiation in Eq. (23) requires an additional boundary 
condition. The most natural condition, which we should 
impose on the solution to Eq. (23), is the requirement of 
boundness at infinity. 
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Equation (19) or, equivalently, Eqs. (23), (24) reduce 
the problem of approximations for U xc to the construc- 
tion of approximations for the xc stress tensor P*£. Since 
the stress density has a clear physical and microscopic 
meaning there is a hope that the later problem is more 
tractable. 



III. STATIC LDA VS. TIME-DEPENDENT LDA 

Let us first derive the standard static LDA from the 
force definition of U xc , Eq. (19). Formally the static 
?7 xc DA (x) is the solution to Eq. (19) in the lowest order of 
the gradient expansion. This solution is obtained by in- 
serting for a homogeneous system of the density n(x) 
into the right hand side of Eq. (19). In the homogeneous 
system the stress tensors P^ u and T^ v are diagonal 



P^[n] = 5^P(n), 



S^Po(n), 



where P and Po are the pressure of the interacting system 
and the pressure of an ideal gas respectively. Therefore 
to the lowest order in the density gradients we get 



CW(x) 



S P 

u flf ± xc 



(n(x)), 



(25) 



where P xc = P — Po is the xc pressure of the homogeneous 
system. Substituting Eq. (25) into Eq. (19), and using 
the common thermodynamic relations, dP — nd\x, \x = 
dF/dn, we find the following result for the xc potential 



f/ xc [n](x) = C/ x L c DA (n(x)), U^ A (n) = ^ 



(26) 



Here P xc is the xc free energy of the homogeneous system. 
The result of Eq. (26) recovers the standard static LDA 5 . 

Physically the above derivation of the static LDA 
means the following. If the density distribution is a semi- 
classically slow function in space, every small volume el- 
ement can be formally considered as an independent ho- 
mogeneous many-body system. The density in this ho- 
mogeneous system equals to the density at the location 
of the element. By solving the homogeneous many-body 
problem we find the stress tensor, which, after the sub- 
stitution into Eq. (19) provides us with C/ XC DA . 

The situation in the time-dependent theory is much 
more complicated. Even if at any instant t the density 
distribution n(x, t) is a slow function in space, a small 
volume element, located at some point x, can not be 
considered as system that is independent of surrounding 
space. For a nonadiabatic dynamics, particles, arriving 
at the point x from other regions, bring an information 
about other spatial points. This is the physical reason 
for the well known nonlocality, immanent to any nona- 
diabatic TDDFT 14 ' 15 . It is straightforward to demon- 
strate the failure of any plain attempt to extend the above 
derivation of the static LDA to the time-dependent case. 
Indeed, a homogeneous many-body problem, which we 
would get by formally separating a small volume clement, 



corresponds to an infinite system with strongly noncon- 
served number of particles. Apparently this problem is 
meaningless. 

In the rest of this paper we show that the nonlocality 
problem in TDDFT is resolved by changing a "point of 
view" on the noncquilibrium many-body system. Any 
flow in the system can be considered as a collection of 
small fluid elements moving along their own trajectories. 
It is possible to divide the system into elements in such a 
way that the number of particles in every element will be 
conserved. Indeed, by the proper deformation and rota- 
tion of a fluid clement one can always adjust its shape to 
the motion of particles and thus prevent the flow through 
its surface. Let us attach a reference frame to one of 
those moving elements. The motion of the origin of this 
frame compensates the translational motion of the fluid 
element. By properly changing scales and directions of 
coordinate axes we can also compensate both the defor- 
mations and the rotation. This means that an observer 
in the new frame will see no currents in the system, and 
a stationary density distribution. Thus from the point of 
view of the co-moving observer the noncquilibrium sys- 
tem looks very similar to the equilibrium one, as it is 
seen by a stationary observer in the laboratory reference 
frame. This similarity is of course not complete since par- 
ticles in the described co-moving frame should experience 
inertia forces. However the inertia forces are determined 
only by local geometric characteristics of the frame. The 
locality of inertia forces and the stationarity of the den- 
sity allow us to consider a small volume element in the 
co-moving frame as an independent many-body system. 
Therefore we can extend the derivation of the static LDA 
to the time-dependent case. 

The description of a flow in terms of trajecto- 
ries of small liquid element is the main idea behind 
the Lagrangian formulation of the classical continuum 
mechanics 24 . One can show that the transformation to 
the Lagrangian coordinates exactly corresponds to the 
transformation to the co-moving reference frame. In the 
next section we apply the general description of quantum 
many-body dynamics in the Lagrangian frame 12 to the 
corresponding reformulation of TDDFT. 



IV. MANY-BODY THEORY IN THE 
LAGRANGIAN FRAME AND GEOMETRIC 
FORMULATION OF TDDFT 

A. Many-body problem in the co-moving frame 

First we briefly review the key results of the many 
body-theory in the Lagrangian frame (all details and 
derivations can be found in I) . The co- moving Lagrangian 
reference frame is defined as follows. Let v(x, t) = 
j(x, t)/n(x, t) be the velocity of the flow. By solving the 
following initial value problem 



0x(£,f) 
dt 



= v(x(£,t),i), x(£,0)=£ 



(27) 
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we find the function x(£,i), which describes the trajec- 
tory of a fluid clement. The initial point, £, of the tra- 
jectory can be used as a unique label of the element. 
This initial position of an infinitesimal fluid element is 
called the Lagrangian coordinate. The transformation 
from the original x-space to the £-space of initial posi- 
tions is the transformation from the Eulerian to the La- 
grangian description of a fluid 24 . On the other hand, the 
equation x = x(£,t), which maps x to £, exactly corre- 
sponds to the transformation to the frame, attached to a 
fluid element. One of the most important characteristics 
of the Lagrangian frame is Green's deformation tensor 24 , 
SV(£,*), 



dx a dx a 



dx^dx^ 



(28) 



Tensor g^ v plays a role of metric in the Lagrangian £- 
space (we will use the notation g for the determinant of 
9nv)- It has been shown in I that the field operators, 
i>(£,i), in the Lagrangian frame are related to the field 
operators, ip(x, t), in the laboratory frame as follows 

^(£,t) = s^(x(€,t),t). 

Apparently the operators tp(£,t) satisfy the common 
equal-time commutation relations, which is guaranteed 
by the factor g* in their definition. The current oper- 
ator, j M (£, t), and the density operator, n(£,t), in the 
Lagrangian frame are defined by the following expres- 
sions 



—i ( ~i dip 
2^ I W ~ ~d¥ 



(29) 
(30) 



where v v = 



is the covariant component of the ve- 



locity vector v, transformed to the new frame 



"dx" 



v"(x(£,t),t) 



The Heisenberg equation of motion for the density oper- 
ator of Eq. (29) takes a form of the operator continuity 
equation 



dn djv 



(31) 



On the level of expectation values Eq. (31) is trivially 
satisfied. One can check by the explicit calculations that 
the expectation value of the current operator, Eq. (30), is 
zero, while the expectation value of the density operator, 
Eq. (29), is time-independent 



i"«,*) = 0"' i (€,t)> = o. 

n(Z,t) = $(£,t)) =n(€,0) 



(32) 
(33) 



where n (x) is the initial density distribution. Equations 
(32) and (33) are in complete agreement with the quali- 
tative discussion in the previous section. 

According to the results of the paper I, the vector 
m *v(£>*) pl a y s a ro l e of an effective vector potential in 
the equation of motion for the field operator, ip. In gen- 
eral the velocity vector v^(^, t) has both potential (longi- 
tudinal) and rotational (transverse) parts. The potential 
part of a vector potential can be always removed from 
the kinetic energy operator by the gauge transformation. 
Therefore it is convenient to separate explicitly the po- 
tential part, VLfi — d^ip, of the vector 



dip „ 



(34) 



where vt\i is the transverse part of v^. Performing the 
gauge transformation ip — e imv ip' in the equation of mo- 
tion for ip (see Eq. (34) in I) we obtain the following 
equation of motion for the transformed operator ip' 



m- 



(35) 



r 



where An(£,t) — n(£,t) — n(£,t) (the Hartree term is 
included in U = U cxt + Uh). Other notations in Eq. (35) 
are the same as in I: 

d 

is the operator of kinematic momentum, and l^ ^> is 



erywhere rising and lovering of tensor indexes are per- 
formed according to the standard rules, i.e. = g^A" , 
etc.). The deformation tensor and the velocity vector in 
Eq. (35) describe generalized inertial forces in the local 
noninertial reference frame. Tensor g M „ in the kinetic 
energy term produces the "geodesic" force. This inertia 
force is responsible for the motion of a free particle along 



the length of geodesic connecting points £ and (ev- the geodesic in £-space. The velocity «„, which acts as a 
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vector potential in Eq. (35), produces the Coriolis force 
(an effective Lorentz force) and the linear acceleration 
force (an effective electric field). The last term in the 
brackets in Eq. (35) is responsible for the inertia force 
that is related to the kinetic energy of the frame (an 
analog of the centrifugal force). 

Equation (35) is the equation of motion in a reference 
frame moving with some velocity v. In fact, the form of 
Eq. (35) is covariant under an arbitrary transformation 
of coordinates, which is generated by a continuous vector 
valued function v(x, t) via Eq. (27). To specify a par- 
ticular reference frame we need to impose an additional 
"gauge" condition. The gauge condition assigns a par- 
ticular value to the generating function v(x, t). There 
are a few formal possibilities to specify the co-moving 
Lagrangian frame (see I). For example, since the expec- 
tation value of the current operator in the Lagrangian 
frame should be zero, we can impose the condition of 
Eq. (32) on the solutions to the equation of motion, 
Eq. (35). In the present paper we prefer to use another 
gauge fixing condition. Namely we require that the so- 
lution to the equation of motion, Eq. (35), should be 
consistent with the equation of the local force balance in 
the Lagrangian frame 



dv> 



7> 



d ( dtp v v v l 



Vff pv _ r, 

no 

(36) 



where P" is the covariant divergence 34,35 of the stress 



tensor 



pv 



V9 dp 2 3£ 



(37) 



The force balance equation of Eq. (36) corresponds to 
the local momentum conservation law, Eq. (11), trans- 
formed to the Lagrangian frame. The stress tensor in 
the Lagrangian frame, P^, which enters Eq. (36), is a 
linear functional of the one particle density matrix pi, 
and of the pair correlation function G^: 



[IV 



P»v[puG 2 ]{£,t) 



(38) 



The explicit microscopic form of the functional 
P/j,v[pi,G2] is presented in Appendix A (sec Eqs. (Al), 
(A4) and (A5)). 

Equation (36) has precisely the same physical signif- 
icance as the static force balance equation of Eq. (16). 
It shows that the inertia forces exactly compensate the 
external force, -g^U , and the force of internal stress, 

^■P^. v . The result of this compensation is the absence 
of the current, and the stationary density distribution in 
the Lagrangian frame. 

Equations (35) and (36) constitute the full set of equa- 
tions of quantum many-body theory in the Lagrangian 
frame. 



B. TDDFT in the Lagrangian frame. Stress tensor 
as a universal functional of the dynamic deformation 

Now we are ready to the discussion of TDDFT. The 
complete description of many-body dynamics in the La- 
grangian frame corresponds to the solution of the equa- 
tion of motion, Eq. (35), supplemented by the frame fix- 
ing condition of Eq. (36). Let us note that both Eq. (35) 
and Eq. (36) contain the same effective potential (the 
term in the brackets in Eq. (35) and Eq. (36)). Using this 
simple property we can formulate the following two-step 
procedure for solving the system of Eqs. (35), (36). On 
the first step we solve a universal nonlinear many-body 
problem of the form 



dt 



2m ! 



+ 



J deVk€')An(0 + *W€,*) 



V/(|)(39) 



where the effective potential U s - C (£, t) is the solution to 
the following self-consistency equation 



dv- 



m - 



7> 



dt 



(40) 



The initial conditions for Eqs. (39) and (40) are the same 
as in the original physical many body-problem, Eqs. (35), 
(36). The special form of the self-consistency equation, 
Eq. (40), ensures the stationarity of the particles density 
and zero current density. Indeed, using Eqs. (39) and 
(40) we obtain the following equations of motion for the 
density n and for the current j 



dn dj>* 
h — = 

dt dp 



djp_ 
dt 



+ F lw j v = 



(41) 
(42) 



where F M „ is the skew-symmetric vorticity tensor, which 
plays a role of an effective magnetic field 



- 1 hi ; — 



dv Tv dv u 



dp dp dp 



dv v 

dp' 



(43) 



The cancellation of the "external" force VJ7 S _ C , and the 
inertial and the stress forces in Eq. (42) is a consequence 
of the self-consistency equation, Eq. (40). By solving 
Eqs. (41), (42) with the initial conditions j(£,0) = 0, 
and n(£, 0) = no(£), we indeed confirm that for alii > 

j(£,i) = 0and n(€,t)=no(€)- 

The self-consistent nonlinear problem of Eqs. (39), (40) 
is universal in that sense that no external potential enters 
the equations. The only "external" variables in Eqs. (39) 
and (40) are the deformation tensor g^ v (^, t) (an effective 
metric), and the transverse part of the velocity, vt(i(£, t), 
(an effective vector potential). The vector vt^ is uniquely 
determined by the skew-symmetric vorticity tensor 
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(an effective magnetic field), Eq. (43). Therefore by solv- 
ing the nonlinear problem of Eqs. (39), (40) we find the 
many-body density matrix as a functional of two-basic 
variables, g^($,t) and i*^^, f) 36 . Inserting this den- 
sity matrix into the microscopic expression for the stress 
tensor, Eq. (38), we obtain the universal functional 



P 



(44) 



Calculation of the stress tensor functional, Eq. (44) com- 
pletes the first step in the solution of the original many- 
body problem. 

The symmetric Green's deformation tensor, g^, and 
the skew-symmetric vorticity tensor, F M „, completely 
characterize the deformed state of a fluid in the La- 
grangian description. Therefore Eq. (44) can be inter- 
preted as the exact noncquilibrium "equation of state" 
that relates the stress tensor to the dynamic deformation 
in the system. Since g^ v and F^ are the functionals of 
velocity, the stress tensor of Eq. (44) is also a functional of 
velocity in agreement with Runge-Gross theorem. How- 
ever, the present interpretation of P^ v as a deformation 
functional looks more natural physically. 

Substituting the "equation of state", Eq. (44), into the 
force balance equation of Eq. (36) we get the exact quan- 
tum "Navier-Stokes" equation in the Lagrangian formu- 
lation. The full set of the exact hydrodynamics equations 
includes Eq. (36) and the trajectory equation, Eq. (27). 
The solution of the system of Eqs. (27), (36) corresponds 
to the second step in the solution of the original many- 
body problem. On this step we determine the evolution 
of velocity for a given external potential. Equations (27) 
and (36) with the stress tensor of Eq. (44) correspond 
to the exact TDDFT hydrodynamics in the Lagrangian 
formulation of continuum mechanics. 

The KS formulation of TDDFT requires a knowledge 
of the xc potential U xc . In Sec. II we have shown that U xc 
is related to the xc stress tensor = P M „ — T^ v , where 

T^ v is the stress tensor for the noninteracting KS sys- 
tem. Obviously, the KS stress tensor can be found from 
the solution of a nonlinear noninteracting problem that 
corresponds to Eqs. (39), (40) with w{l^>) = 0. Hence 
by solving Eqs. (39), (40) with and without interaction 
we compute P^ v and T^ v respectively. The difference 
of these tensors gives us the xc stress tensor in the La- 
grangian frame as a functional of g^ v and F^ v 



(45) 



Transforming the xc stress tensor of Eq. (45) back to 
the laboratory frame, and substituting the result into 
Eqs. (23), (24), we obtain the equation for the xc po- 
tential C/ xc (x, t). Another possibility is to determine xc 
potential, i7 xc (£,i), directly in the Lagrangian frame by 
solving the following equation 



(46) 



The transformation of U xc (£,t) to the laboratory frame 
corresponds to the following replacement £ — > £(x, t), i.e. 



U xc (x,t) = U xc [ P M „](|(x, t),t) 



(47) 



where U xc is the solution to Eq. (46). 

Let us note that the problem of calculation of the equi- 
librium stress tensor functional, P^u, in the static DFT 
can be formulated in exactly the same fashion. To cal- 
culate P^v [n] (x) we need to find the equilibrium solution 
to the following universal nonlinear many-body problem 

d V 2 

= -^( X ) + ^s-c(x)V(x) 

ot 2m 



+ 



d_ 



J dxX|x-x'|)An(x')V(x) (48) 

(49) 



I d 



— £/ s _ c (x) = -—p^[p u G 2 ], 



where P^ [pi , G 2 ] = T^[pi] + W^[G 2 ] is defined after 
Eqs. (14), (15). For a given density n(x) the equilibrium 
solution to Eqs. (48), (49) defines the stress tensor P^ v 
as a universal functional of n. 

Stationarity of the density in the Lagrangian frame 
makes the dynamic universal problem of Eqs. (39), (40) 
to a large extent similar to the equilibrium universal 
problem of Eqs. (48), (49). In Sec. V we use this sim- 
ilarity to derive a local nonadiabatic approximation in 
TDDFT. 



V. TIME-DEPENDENT LOCAL 
DEFORMATION APPROXIMATION 

In the previous section we have shown that the calcu- 
lation of xc stress tensor, which defines the xc potential, 
reduces to the solution of the nonlinear universal many- 
body problem, Eqs. (39), (40). Obviously, it is not possi- 
ble to solve this problem exactly. However one can try to 
find an approximate solution by a perturbative expansion 
in terms of some small parameter. Below we construct a 
local approximation that corresponds to the lowest order 
in the gradients of basic variables (i.e. the density in the 
static DFT, and the deformation tensor in TDDFT). 



A. General formulation of a nonadiabatic local 
approximation 

1. Preliminaries: Derivation of the static LDA 

To illustrate the general procedure we start again with 
the familiar case of the equilibrium theory. The prob- 
lem is to find the equilibrium solution to Eqs. (48), (49), 
assuming that gradients of the density n are vanishingly 
small. In the limit Vn — > the spatial derivatives of the 
stress tensor also vanish. Hence to the lowest order in 
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the density gradients the solution to the self-consistency 
Eq. (49) takes a trivial form, f7 s - c (x) = G, where G is a 
constant. Therefore the many-body equation of motion, 
Eq. (48), simplifies as follows 



.dip(x) 



v 2 



dx'iy(|x-x'|)ri(x / )V'(x). (50) 



Thus the nonlinear problem of Eqs. (48), (49) reduces to 
the usual linear many-body problem for a homogeneous 
equilibrium system with a given density n. Substituting 
the equilibrium solution to Eq. (50) into Eqs. (13)— (15) 
we compute the stress tensor, P^ v {n) = ^ jUJ/ P(n), where 
P is the pressure of the homogeneous system: 



P(n) 



d 



kin 



2d J 



' 9x 



G c 2 q (|x|)dx. 



(51) 



Here d is the number of spatial dimensions, Pkin is the 
kinetic energy per unit volume, and G^dxl) is the pair 
correlation function of the equilibrium homogeneous sys- 
tem. Similarly by solving the homogeneous nonintcr- 
acting problem we find the KS stress tensor T£,(n) = 



^g-Ekin' where P^n is the kinetic energy of an ideal 
Fermi gas. Substituting P jU1/ (n(x)) and T^(n(x)) into 
Eqs. (20) and (19) we recover the common static LDA 
(see Sec. III). 



,(o) 



2. Basic equations of TDLDA: The homogeneous 
many-body problem 

The above procedure allows for a straightforward ex- 
tension to the time-dependent problem. Let us assume 
that the characteristic length scale, L, of the deformation 
inhomogeneity goes to infinity. In this limit the vector 
7^Pp., y in the right hand side in Eq. (40) vanishes. There- 
fore to the lowest order in 1/L — > the self-consistent so- 
lution to Eq. (40) takes the form: J7 s _ c (x, i) = C(t) and 
vt = 0. Substituting this solution into Eq. (39) we get 
the equation of motion for f/'-operator 



. om = g ^) 9 2 m 

1 dt 2m d^df 



dt'w(U-t'\\)WMt) 

(52) 



where — = is the length of geodesic in a homo- 
geneously deformed Lagrangian space (see Appendix B): 



ll€ - a = V-9^W(e - - €'")■ (53) 

Equation (52) corresponds to a homogeneous many- 
body system. It is more natural to reformulate this ho- 
mogeneous problem using the momentum representation 
for field operators 



(54) 



The equation of motion for annihilation operator, ak, 
takes the form 

.da k »„,.nMv~ , \- w(||q||) t 



dt y w 2m 



° k + E — 7^r a P a p+q ak -q ( 55 ) 

D.a V9\P) 



where w{q) is the Fourier component of the interaction 
potential, and 



(56) 



is the norm of the wave vector in the deformed momen- 
tum space. Equation (55), corresponds to the following 
Hamiltonian 



H 



9 ^m"< ak 



1 



2^9 



k,q 



(57) 

where n q = J2 P a P a p+q is the density operator in the 
momentum representation. In Appendix B we show that 
to the lowest order in the deformation gradients the mi- 
croscopic expression for the stress tensor, P M „, simplify 
as follows 



Liu - 



1 



E 

k 



-E 



/(k) 



w 



|k||)+^(||k|| 



G 2 (k), (58) 



where f(k) = (a k &k) is the Wigner function, C?2(k) is 
the Fourier component of the pair correlation function, 
and w'(x) = dw(x)/dx. 

It should be mentioned that Eq. (58) can be derived di- 
rectly from the "geometric" definition of the stress tensor 
(see I and Ref. 37) 

Indeed, using the relation Sg = —gg^vSg* 1 " , and com- 
puting the derivative of the Hamiltonian, Eq. (57), with 
respect to g^ u , we immediately recover Eq. (58). 

The Hamiltonian of Eq. (57) determines the homoge- 
neous problem which we need to solve for the derivation 
of a local approximation in TDDFT. This problem corre- 
sponds to a system of particles in a small volume located 
at the point £ of Lagrangian space. The density of parti- 
cles is time-independent and equals to the initial density, 
n o(0> a t that point (obviously, the operator of the num- 
ber of particles commutes with H) . The behavior of the 
system is governed by the local value of the deforma- 
tion tensor, g^i^, t). By solving the equations of motion 
we find the Wigner function, /(k, t), and the pair cor- 
relation function, (j?2(k, t). Substitution of /(k, t) and 
G 2 (k, t) into Eq. (58) yields the stress tensor functional 

P^Ag^i^, t), n o(0]- By the repetition of the above pro- 
cedure for the noninteracting system (Eqs. (57), (58) with 
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w = 0) we find the KS stress tensor, T^lg^fe, t), no(£)], 
and, finally, the xc stress tensor 



(59) 



Substituting P*£ of Eq. (59) into Eq. (46) we determine 
the corresponding xc potential in the Lagrangian frame. 

The xc stress tensor Eq. (59), is a local in space 
functional of the deformation tensor (it should be noted 
that in general this functional is nonlocal in time). In 
what follows the approximation of Eq. (59) will be re- 
ferred to as a Time-Dependent Local Deformation Ap- 
proximation (TDLDA). 

The construction of TDLDA reduces to the solution 
of the homogeneous many-body problem. In this re- 
spect the situation is similar to the static case. How- 
ever, the homogeneous time-dependent problem, defined 
by the Hamiltonian of Eq. (57), is still too complicated 
to be solved exactly. Indeed the operator equation of 
motion, Eq. (55), generates an infinite set of coupled 
evolution equations (BBGKY hierarchy 38 ) for correlation 
functions. The first equation of this hierarchy is the equa- 
tion of motion for the Wigner function 



.dm 
dt 



V9 



a l ( a k fl k-q 



Z lc+q 



An equation for the four-fermion correlator, entering the 
right hand side in Eq. (60), couples to the six-fermion cor- 
relation functions, etc. However the homogeneity of the 
problem and a very specific form of the "driving force" in 
the equations of motion allow us to construct reasonable 
approximate xc functional (see next subsections). 



3. Stress tensor of the noninteracting KS system 

A necessary step in the derivation of TDLDA is to com- 
pute the stress tensor, T^, in the noninteracting system. 
This problem can be solved exactly. In the noninteracting 
case (w = 0) Eqs. (58) and (60) reduce to the following 
simple form 



f s = — V 



^/(k,t), 

m 



d_ 

dt 



/(k,i) = 0. 



(61) 
(62) 



Equation (62) shows that the distribution function of 
noninteracting particles in the Lagrangian frame is time- 
independent. Let us assume for definiteness that the sys- 
tem evolves from the equilibrium state, i. e. /(k, 0) = 
where n£" is the Fermi function. In this case the solution 
to Eq. (62) takes the form 



/(k,t)=/(k,0) 



(63) 



Substituting Eq. (63) into Eq. (61) we get the kinetic 
stress tensor of the KS system in the Lagrangian frame 



5p V 



V9W) 



(64) 



where the function Pq{u) = | E^L( n ) IS * ne equilibrium 
kinetic pressure of a noninteracting homogeneous Fermi 
gas. 

For the practical calculation of xc potential in TDLDA 
(see Eqs. (19), (20) ) we need the stress tensor, T^(x, t), 
in the laboratory frame. Application of the common ten- 
sor transformation rules 34 , 

P^t) = ^^P a ^,t),t), (65) 
to the stress tensor of Eq. (64) yields the result 



TL (x, t) = (x, t) V^KOPo (no (£ (x, t) ) , (66) 



were g^{x.,t) is the Cauchy's deformation tensor 24 



6VCM) 



dx^ dx u 



(67) 



The determinant, g(x,t), of Cauchy's deformation ten- 
sor, Eq. (67), is related to the determinant, g(£,t), of 
Green's deformation tensor, Eq. (28), as follows 



g(x,t)=g- 1 (Z(x,t),t). 



(68) 



Equation (66) clearly demonstrates an extreme nonlocal- 
ity which is related to the memory effects. The stress 
tensor, T^(x, t), at a given point x depends on the ini- 
tial density at the point £(x, t) that is the initial position 
of a fluid element presently at x. Let us show that this 
dependence on the delayed position can be represented 
in a local form. By definition of the Lagrangian coordi- 
nate, the density, n(x, t), in the laboratory frame can be 
expressed in terms of the initial density distribution (the 
density in the Lagrangian frame): 



n(x, t) 



n (g(x, t)) 



(69) 



Using the relation of Eq. (68) we can represent the non- 
local quantity n (£(x, t)) in the following form 



n (£(x,i)) 



n(x, t) 



(70) 



Substituting Eq. (70) into Eq. (66) we obtain a com- 
pletely local representation for the KS kinetic stress ten- 
sor 

The nonlocality of the stress tensor in the form of Eq. (66) 
is now hidden in the space-time dependence of the func- 
tion g(x, t). 
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B. Exchange-only TDLDA 

The most difficult part in the derivation of an explicit 
TDLDA is the solution of the interacting problem defined 
by the Hamiltonian of Eq. (57). In this subsection we find 
the exact solution of this problem in the exchange ap- 
proximation, which provides us with the x-only TDLDA. 
In the x-only case the pair correlation function G2 (k, t) 
is completely determined by the one particle distribution 
function /(k, t) 



G 2 (k,i) = -^/(k + p,i)/(p,i) 



(72) 



Performing the mean field decoupling of the four-fermion 
terms in Eq. (60) we find that the right hand side in this 
equation vanishes. Therefore the equation of motion for 
the function /(k, t) takes the form 



d_ 

dt 



/(k,i) = 0. 



(73) 



Equation (73) coincides with the corresponding equation 
of motion for the noninteracting system, Eq. (62). Hence 
in the x-only approximation both the Wigncr function 
and the pair correlation function in the Lagrangian frame 
preserve their initial form 

7(M)=n£, (74) 
G 2 (k,t)=GS(no;fc) = -£n£ +p n£. (75) 



where we introduced a shortcut notation w p = w(j>). 
Equations (78), (24) and (23) uniquely determine the lo- 
cal potential U x {x.,i) in x-only TDLDA. Apparently the 
exchange potential U K (x, t) is a local (both in space and 
in time) functional of the density n(x, t) and Cauchy's 
deformation tensor (^(x, t). In the equilibrium system 
(ffjui/ = ^ixv) the potential, defined by Eqs. (78), (24) and 
(23), reduces to that in the common static local exchange 
approximation. 



C. Linear response TDLDA. 

In the linear response regime the deformation tensor, 
g^ v , slightly deviates from the Kronecker symbol 

<r«,t) « <w &<r = -|p-fp (79) 

where u = x — £ is the displacement vector, which is 
assumed to be small. In the linearized theory the trajec- 
tory equation of Eq. (27) reduces to the common linear 
relation of the velocity to the displacement 



dt 



v(€,t). 



(80) 



Substituting g» v of Eq. (79) into Eq. (57) and keeping 
only linear in <5g M1/ terms, we obtain the following lin- 
earized Hamiltonian 



H = H + P^W, 



(81) 



Here k = |k| = y/k^k^ is the usual modulus of k, and 
Gf (n; k) is the exchange pair correlation function in the 
equilibrium Fermi gas of the density n. Substituting 
Eqs. (74) and (75) into Eq. (58) we obtain the follow- 
ing stress tensor in the interacting system 

P M „ = -%=P (n (O) + (no(€),M€.*))- (76) 

The first term in the right hand side in Eq. (76) is the 
kinetic stress tensor of the noninteracting system, while 
the second term, P* v , corresponds to the exchange con- 
tribution to the local stress density 



^ 2g ^ 



^^(||k||)+ 3 ^(||k|| 



Gl{n -k) 



(77) 

Using the transformation rules of Eq. (65) we get the 
following expression for the exchange stress tensor in the 
laboratory frame 



w' p + S^w p 



xG% ( — ; J , (78) 



where H is the standard Hamiltonian for the homoge- 
neous system, and P M „ is the stress tensor operator 



p» v = 2^ [~^~ a k a * + 2 \ l~ Wk f ii/Wk J nkn - k 

(82) 

First we need to compute the stress tensor P^ v in the 
Lagrangian frame, Eq. (58) . In the linear regime Eq. (58) 
takes the form 



P»» = S^P(n ) + Q^H^H 



(83) 



where the linear response kernel, Q^a/si^), can be rep- 
resented as follows 



(84) 



The first frequency independent term, Q^ va p , in Eq. (84) 
comes from the explicit local in time dependence of the 
integrals in Eq. (58) on the deformation tensor g^it)- 
Namely, the fourth-rank tensor Q^ uaf3 is defined by the 
following derivative 



Q 



l±va.(3 



dP^[t\G e *] 



8g 



a/3 



(85) 
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where P M1 ,[/ cq , G^] is the stress tensor of Eq. (58), calcu- 
lated with the equilibrium Wigncr function, / cq (k), and 
the equilibrium pair correlation function, G^k). The 
perturbation (the second term) in the linearized Hamilto- 
nian of Eq. (81) induces deviations of the Wigner function 
and the pair correlation function from their equilibrium 
values. These deviations are responsible for the second, 
nonlocal in time term in Eq. (84) 



kjjty Sf(k,t) 
2m 6g a P(t') 



2 ^ \k~ Wk ^ Wk ) Jg^W)' (86) 

Comparing Eqs. (81), (82) and (86) we find that the 
dynamic kernel AQ^ap can be related to the following 
stress autocorrelation function 



Jo 



dte wt ([P^(t),P a p(0)])- (87) 



Substituting the stress tensor in the Lagrangian frame, 
Eq. (83), into the transformation formula of Eq. (65) we 
compute the stress tensor in the laboratory frame 



P M „(x,u;) = ^P(n (x)) + SP^fcu), 



(88) 



where 



dP 

'dn 



5P^ = - V^uVno + Q^ a p(u)g ap {u). (89) 



The kernel Q^ a p{^) in the laboratory frame takes the 
form, which is similar to that of Eq. (84) 



(90) 



The frequency independent term, Q c ^ va p 1 in Eq. (90) is 
related to the quantity Q^, Q(3 of Eq. (85) as follows 



CO 



(91) 



It is worth mentioning that the first term in the right 
hand side in Eq. (89) guarantees that the harmonic po- 
tential theorem 16 is satisfied. In fact, only this term sur- 
vives for the rigid motion of the system. Within the 
present formalism this term comes from the expansion 
of the argument of P^ in the transformation rule of 
Eq. (65). The correction to the kernel in the labora- 
tory frame (the first term in Eq. (91)) corresponds to the 
expansion of the tensor prefactor in Eq. (65). 

By symmetry the fourth-rank tensor Q^ap is uniquely 
representable in the form 

Q l ivaP = (-^ + ^ j <W<W + 7^(<W^W3 + tW<W); (92) 



The scalar coefficients, K(u>) and fi(u>), in Eq. (92) are 
related to the tensor Q^ a p contracted over different cou- 
ples of indexes 



K(uj) = -pQ aaf3l 3(Lj), 



d 2 + d-2 



Q a pap{u) - -Qaapp(u>) 



(93) 
.(94) 



Substitution of Eq. (92) into Eq. (89) yields the following 
result for the linear correction to the stress tensor in the 
laboratory frame 



SP^ = - S^ v - — uVn 



dP ^ 
dn 
1 



.6 ia ,K^6g aa + n[6g'»'--%-6g ac 



Spy 
d 



(95) 



The stress tensor SP^ of Eq. (95) has a clear visco-elastic 
form, where K(ui) and fi(u>) arc the bulk modulus and the 
shear modulus respectively. The xc stress tensor, 5P™, is 
the difference of the expressions given by Eq. (95) for the 
interacting and the noninteracting systems. Apparently 
5Pj£ takes a form of Eq. (95) with P, K and ft being 
replaced by P xc , K xc and ^ xc respectively, where 



-P (0) 
d ^ 



K XC = K-K ( 



()■ 



Mxc = M — Mo i 



(96) 
(97) 



are the xc pressure and the xc visco-elastic moduli. In 
Eqs. (96), (97) p£5n, Ko, and fiQ correspond to the ki- 
netic energy, the bulk modulus, and the shear modulus 
of an ideal Fermi gas. Therefore in the linear response 
regime our TDLDA naturally reduces to the Vignale- 
Kohn approximation 17 in the visco-elastic formulation of 
Ref. 18. 

An explicit microscopic representation for the bulk, K, 
and the shear, fi, moduli can be found using Eqs. (94), 
(93), (91), (90), (87) and (85). Both K and fi take the 
following general form 



K(w) = K°° + AK(u), 
fi(uj) = n°° + Afi(cj). 



(98) 
(99) 



The first terms, K 00 and fi°°, in the right hand sides in 
Eqs. (98) and (99) are obtained by the substitution of 
Q~ Q/3 , Eq. (91), into Eqs. (93) and (94). Performing 
straightforward calculations for the interacting and the 
noninteracting systems we arrive at the following results 
for the high frequency parts of the xc elastic moduli 



K c 



d + 2 



k 2 w% + (3d + l)kw' k + 2d 2 w k 



2d{d + 2) 



G c 2 \k) , (100) 



Mxc — ^kin 



(d+l)kw' k 



2d(d + 2) 



±Gl\k), 



(101) 
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where E\ 



kj,, — -Ekin — is the xc kinetic energy of 
the equilibrium system. In the special case of Coulomb 
interaction the momentum integrals in Eqs. (100), (101) 
can be expressed in terms of the potential energy per unit 
volume, 



In d dimensions the Coulomb potential is proportional to 
l/fc^" 1 . Therefore we get the following identities for the 
derivatives, which enter Eqs. (100), (101) 

kw' k = —(d — l)wk and k 2 w k = d(d — l)wk- 

These identities allow us to simplify Eqs. (100), (101) as 
follows 



2(d + 2) 



fxc — ^ kin 



E kin + 

d-l 
d{d + 2) 



d+l 



pot, 



E, 



pot • 



(102) 
(103) 



The high frequency forms of Eqs. (102), (103) are well 
known in the literature 39 ' 40,41 . Commonly they are de- 
rived using the "third moment sum rule". Within our 
formalism the expressions of Eqs. (100) and (101) for 
and /i^ come about almost trivially from the explicit lo- 
cal in time dependence of the stress tensor, Eq. (58), on 
the deformation tensor. 

To represent the frequency dependent parts of vis- 
coelastic moduli in the most convenient form we decom- 
pose the stress tensor operator P M1/ , Eq. (82), into a scalar 
and a traceless parts 

p iu , = 5 iu ,p + k iu ,, ( 104 ) 

where P = iTr_P M!/ is the pressure operator 

P = 2 ^ |_^"k«k + ^{kw' k + dw fe )n k n_k (105) 

k 

and tt^ v is the operator of the traceless part of the stress 
tensor (Tr7r M „ = 0) 



kw' k + dw k 
2d 



n k n_ k 



k \rgr--d ){m a ^ + 

(106) 

Substituting AQ^^w) of Eq. (87) into Eqs. (93), (94), 
and using Eqs. (104)-(106) we find that AK xc (lo) equals 
to the pressure autocorrelation function 



00 

Ah\,U) = 2i J,'- ( /'(/). /'(())] (107) 



whereas the dynamic part of the shear modulus, 
A/u xc (w), is proportional to the autocorrelation function 
of the traceless tensor operator tt^ v 

00 

A^xcH = d2 2 / e^([9, u (t)^ u (0)})dt (108) 



The coefficient, 2/(d 2 + d-2), in Eq. (108) is exactly the 
inverse number of independent components of a second- 
rank traceless tensor. We would like to outline a very 
natural form of Eqs. (107) and (108), which is in clear 
agreement with the physical significance of the quantities 
K and /1. 

The frequency dependent contributions to the visco- 
elastic moduli are related to the dynamics of the Wigner 
function and the pair correlation function (see Eq. (86)). 
In Sec. VB we have shown that in the exchange approx- 
imation the time-dependent deformation tensor in the 
Hamiltonian of Eq. (57) does not induce any dynamics of 
/(k) and G*2(k). Therefore only correlations are respon- 
sible for nonvanishing Aif xc (cj) and A/U xc (w). The most 
important effect of the dynamic correlations, which is de- 
scribed by Eqs. (107) and (108), is the memory loss due 
to collisions. Since in a zero-temperature Fermi system 
the collisional dissipation should be suppressed, there is a 
hope that the frequency dependent parts of K xc and /x xc 
do not substantially influence the dynamics. [We note 
that this is not in general true for steady state transport 
situations, where the dissipation plays an essential role.] 
Neglecting AX IC (w) and A / u xc (w) we get a purely elastic 
xc stress tensor with the bulk and the shear moduli de- 
fined by Eqs. (100) and (101). Another argument due to 
Conti and Vignale 39 also shows that for an electron gas 
the elastic approximation should work reasonably well. 
Indeed, the dissipation effects are absent in the x-only 
approximation that is valid in the weak coupling (high 
density) regime. In the strong coupling (low density) 
limit electrons tend to form a Wigner crystal - the state 
where the collisional dissipation also vanishes. Therefore 
one naturally expects that at all intermediate densities 
the purely elastic approximation should provide a rea- 
sonable description of the dynamic stress. 



D. Nonlinear elastic TDLDA 

The linear VK approximation with purely elastic 
bulk modulus K^, Eq. (100), and shear modulus fi™, 
Eq. (101), allows for a simple nonlinear extension. In 
this subsection we derive this nonlinear clastic TDLDA, 
and formulate a complete set of self-consistent KS equa- 
tions in a convenient for practical applications form. 



1. Exchange-correlation stress tensor in the elastic TDLDA 

Elastic TDLDA is based on the assumption that both 
the Wigner function and the pair correlation function in 
the Lagrangian frame preserve their initial form. This 
corresponds to the dynamics with extremely pronounced 
memory that is not destroyed by the effects of colli- 
sional relaxation. To get the stress tensor for the system 
evolving from the equilibrium state we have to substi- 
tute J(M) = / cq (n ;fc) and G 2 (k,i) = G c 2 q {n ;k) into 
Eq.(58). As a result the stress tensor in the Lagrangian 



14 



frame takes the form 



+ g^w(\\H) G c 2 q (n ;fc). 



(109) 



Tensor P^^, t) of Eq. (109) locally depends on the den- 
sity, no(£), and Green's deformation tensor, g^ v ($,t), in 
a given point £ of the Lagrangian space. Transforming 
this tensor back to the laboratory frame and subtracting 
the KS stress tensor of Eq. (71) we obtain the following 
result for the xc stress tensor in the physical x-space 



pxc _ ^9^v 

d 



+ 5 llu w p 



(110) 



Equation (110) determines the xc stress tensor as a func- 
tion of the time-dependent density, n(x, t), and Cauchy's 
deformation tensor g^ v (x, t) . Let us remind that the 
memory related nonlocality of P*£, Eq. (110), is hidden 
in the "local" representation of the function no(£(x, i)) 
(see Eq. (70)). The "elastic" xc potential is the solution 
to the Poisson Eq. (23), where the xc "charge density" is 
defined after Eq. (24) and (110). 

In the exchange approximation the stress tensor P*£, 
Eq. (110), reduces to the x-only tensor P* of Eq. (78) 
that is exact in the weak coupling limit. In the linear 
response regime the corrections to the density and to 
the Cauchy's deformation tensor are proportional to the 
displacement vector 



n = n -Vn u, = 5^ - ^ - (111) 



du^ du v 
dx u dx^ 



Linearizing the stress tensor of Eq. (110) and 
using Eq. (Ill) we straightforwardly recover VK 
approximation 17,18 with the elastic moduli K™, 
Eq. (100), and Eq. (101). 



2. Self-consistent Kohn-Sham equations 

Let us formulate the complete set of self-consistent KS 
equations in the elastic TDLDA. The Kohn-Sham formu- 
lation of TDDFT allows to calculate the density n(x, t) 
and the velocity v(x, t) in the interacting iV-particle sys- 
tem using the ideal gas formulas 



N 



n(x,i)=]T|^(x,i)| 2 , 



(112) 



1 N ' 

v ( x ' *) = n E 2^ [^' V ^ ~ ■ (113) 

j'=i 



Single particle orbitals <j>j (x, t) satisfy the time- 
dependent KS equations 



dt 



+ (t^ext + U eS [n,g^ v ])< 



(114) 



where U ext (x, t) is the external potential. For the practi- 
cally important case of 3D system with Coulomb interac- 
tion the effective potential U e s[n, <7 M „](x, t) is the solution 
to the following Poisson equation 



V 2 [/ C ff = An(e 2 n + p xc [n, g^ v ]). 



(115) 



The first term in the brackets in Eq. (115) generates the 
Hartree potential, Un, while the second term is respon- 
sible for the xc potential. The xc "charge density", p xc , 
is the local functional of n and g^ u 



1 d 
4^dx^ 



(116) 



where P™(n, g^ v ) is the function of n(x, t) and (/^(x, t), 
which is defined by Eq. (110). In Appendix C we show 
that for a Coulomb system Eq. (110) simplifies as follows 

9 {117) 

where L^(g a p) is a purely geometric factor that is explic- 
itly defined in Appendix C. Therefore the dependence of 
P*£(n, g^v) on g M „ and on nj \fg is completely factorized, 
which should significantly simplify practical applications. 
The kinetic, E£f n (n), and the potential, E pot (n), ener- 
gies of the homogeneous electron gas can be expressed 
in terms of the xc energy per particle, e xc (n) (see, for 
example, Ref. 39). For d = 3 we get 



£$?(n) = 3n* P£ 
,n3 



E, 



■pot(n) = -3na — 

713 



Hence our nonadiabatic TDLDA requires only a knowl- 
edge of the function e xc (n) for the homogeneous electron 
gas, exactly as the common static LDA does. 

The density n, which enters Eqs. (115)— (117), is related 
to KS orbitals via Eq. (112). The second basic variable, 
Cauchy's deformation tensor g^,, is uniquely determined 
by the velocity v(x, t), Eq. (113). To compute the de- 
formation tensor we need to solve the trajectory equa- 
tion of Eq. (27) and then substitute the solution into the 
definition of g^ u , Eq. (67). It is, however, more conve- 
nient to determine this tensor directly from the solution 
of an equation of motion for g^ u (x,t) 42 . This equation 
of motion can be derived as follows. Let us consider the 
contravariant tensor g^" (the inverse of g^) 



dx^ dx v 



(118) 



Using the trajectory equation of Eq. (27) we can compute 
the time derivative of g^ u , Eq. (118), at constant £ (i.e. 
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within the Lagrangian description) 
'dg» v \ dv» 



dt 



dx a 



9 



dx a 



(119) 



The time derivative of g^ v can be related to the time 
derivative oig >lv = (g^) -1 as follows dtg = —g{dtg~ 1 )g. 
Using this relation and taking into account that 

(dth = (d t )* + vV, 

we get the final equation of motion for Cauchy's defor- 
mation tensor g^(x, t) 



ggy _ a dg^ _ dv a _ dv 

V „ _ „ ..gav 



dt 



dx a dx^- 



dx 



(120) 



Equation (120) should be solved with the initial condition 
<7 M „(x, 0) = 8fj, v , which follows from the initial condition 
for the trajectory equation of Eq. (27). 

The system of Eqs. (112)-(117), (120) constitute the 
complete set of self-consistent KS equations in the nonlin- 
ear elastic TDLDA. In the equilibrium situation {g^ = 
Sf^u) this system reduces to the common static KS equa- 
tion with the LDA xc potential. In the linear regime it 
recovers the results of VK approximation with the elas- 
tic moduli of Eqs. (102), (103). The nonadiabatic mem- 
ory effects are described by Cauchy's deformation tensor, 
which satisfies Eq. (120). It should be noted that from 
the computational point of view the solution of this equa- 
tion do not introduce any addition difficulties. Formally 
Eq. (120) has the same structure as the time-dependent 
KS Eq. (114). Hence Eqs. (114) and (120) can be solved 
simultaneously by the same method. 

Very recently VK approximation has been success- 
fully applied to the description of optical and polariza- 
tion properties of many different systems, such as atoms, 
molecules, semiconductors and polymers 43 ' 44,45 ' 46,47 . 
Since VK approximation is a linearized version of our 
theory, we hope that the general TDLDA also will be- 
come a useful tool for studying nonlinear time-dependent 
phenomena. 



VI. CONCLUSION 

TDDFT extends powerful ideology of the ground state 
DFT to the domain of noncquilibrium phenomena. How- 
ever, in contrast to the static DFT, which is currently a 
common computational tool in many branches of physics, 
its time-dependent counterpart still suffer from a number 
of unresolved problems. One of those problems is a lack 
of well founded basic local approximation that would play 
a role similar to LDA in the static DFT. In this paper we 
have shown that the local approximation in TDDFT can 
be regularly derived, but this derivation requires almost 
complete reconsideration of the theory. We reformulated 
TDDFT from the point of view of a local co-moving ob- 
server. The new formulation of the theory shows that the 



most natural basic variables in TDDFT are the local ge- 
ometric characteristics of the deformations in a quantum 
many-body system. 

Throughout this paper we used the analogy of TDDFT 
to the classical continuum mechanics. The importance of 
the hydrodynamic interpretation, which perfectly fits the 
very idea of DFT, is one of the messages of the present 
work. Using the hydrodynamic formulation of TDDFT 
we were able to relate the xc potentials to the local stress. 
In particular we proved that the exact xc force must have 
a form of a divergence of a second-rank tensor. The well 
known zero force and zero torque sum rules are direct 
consequences of this strong local requirement. The func- 
tional dependence of xc potential on the basic variables 
also acquires a clear physical meaning. It corresponds 
to the stress-deformation relation, which is very natural 
from the point of view of continuum mechanics. If spa- 
tial derivatives of the deformation tensor are small, the 
stress-deformation relation becomes local and therefore 
we get the local approximation for the xc potential in 
TDDFT. It is natural to abbreviate this approximation 
as TDLDA, which means time-dependent local deforma- 
tion approximation. In the linear response regime the 
general stress-deformation relation (TDLDA) reduces to 
the linear Hook's law 30 , which exactly coincides with the 
visco-elastic VK approximation 17,18 . The formal appli- 
cability conditions for the general nonlinear TDLDA are 
the same as for the linear VK approximation. It should 
be also noted that neither of previously proposed nonlin- 
ear phenomenological constructions 18 ' 19,20 is in general 
confirmed by the present regular microscopic considera- 
tion. 

In the last section of this paper we introduced the clas- 
tic TDLDA. In this approximation the xc stress tensor is 
simply a function of the density and of the Cauchy's de- 
formation tensor. For a system with Coulomb interaction 
we presented the xc stress tensor and the xc potential in 
an explicit "ready for implementation" form. We also 
formulated the full set of self-consistent KS equations in 
TDLDA. In the equilibrium state the deformation ten- 
sor is diagonal and TDLDA reduces to the standard 
static LDA, while in the linear response regime it recov- 
ers VK approximation. To conclude we mention that the 
self-consistent equations of Sec. VD2 can be straightfor- 
wardly reformulated in terms of xc vector potential. The 
only difference is that the Poisson equation for U xc should 
be replaced by Eq. (21) which relates A xc to the xc stress 
tensor. This replacement introduces one more evolution 
equation which should be solved simultaneously with the 
KS equation and the equation for Cauchy's deformation 
tensor. 
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APPENDIX A: STRESS TENSORS IN A 
GENERAL NONINERTIAL FRAME 

The microscopic representation for the stress tensor 
P[i.v{£,,i) in a general local noninertial frame has been 
derived in I: 

P^t)=T^,t) + W^,t), (Al) 

where the kinetic stress tensor, T^„(£,t), and the inter- 
action stress tensor, W pl ,(£, t), are obtained by the trans- 

I 



formation of T^x, t), Eq. (14), and W^x, t), Eq. (15), 
to the new frame. Namely, 

^(€,*) = ^^Ta/3(x«,t),t), (A2) 
t) = ^^fw aP (x(t t), t). (A3) 

The result of the transformation, Eqs. (A2), (A3), takes 
the following form 



WV(£,i) = — ^— ^ / cZA / drjdrj 6($- z„,„>(A)) — ^ G 2 (v,V )■ (A5) 

•'0 j c/>"n,ri' 



Here the function z,j jT? /(A) is the geodesic that connects 
points 77 and 77', and l v , v > is the length of this geodesic. 
The curve Zjj.jj' (A) can be found from the solution of the 
geodesic equation (see, for example, Ref. 34) 



^(A) + r^(z)i«(A)^(A) = 0, 



(A6) 



supplemented by the boundary conditions z(0) = 77, 
z(l) = 77'. In Eqs. (A5) and (A6) z = dz/dX, and T» is 
affinc connection: 



1 uu ( d 9u a . dg v p dg af} 



2 y d£ a d£ u 



(A7) 



Equations (A4) and (A5) define tensors T^ v and W M „ 
as functionals of the microscopic state of the system. 
Tensor is a linear functional of the one particle den- 

sity matrix, pi (£,£') = $t(€MO>- Similarly, 

is a linear functional of the pair correlation function 

- (^(€)n(€'W€)> -n(€)«(0- Therefore 
P„„ = + WV = G 2 ](€, t) (A8) 

Equation (A8) is the result, which we need for the dis- 
cussion of TDDFT in Sec. IV. 



APPENDIX B: STRESS TENSORS FOR A 
HOMOGENEOUS DEFORMATION 

For a homogeneous system with <7 M „(£, t) = <7 M „(i) and 
VTfj, = the general expressions, Eqs. (A4) and (A5), for 
the stress tensors simplify as follows. Equation (A4) for 



the kinetic stress tensor of takes the form 

1 / dip' dip dip'' dip 



T 



2m^/g \ d£» d£ v df d^ 

1 3 2 Pi(g) 

rn^/g dfdf ' 



(Bl) 



where pi(£ — £') = (V'^(€)V'(€')) i s the one particle den- 
sity matrix for the homogeneous system. Introducing the 
Wigner function 

/(k) = J e- ifc ^Vi(€)d€, (B2) 

we obtain the following final representation for T M „ 



(B3) 



To calculate the interaction stress tensor, Eq. (A5), we 
need to solve the geodesic equation of Eq. (A7). For a 
homogeneous metrics the solution is a straight line: 



z r,,r,'(A) =T7+ (77'- T7)A. 



(B4) 



The length of the geodesic, l v ,r)>, can be calculated as 
follows 

l n,* = f ^fg^(z)zn(\)z»(\)d\ 

= yjg^rT - v»){v" ~ ff) ■= \\V ~ V'l (B5) 

Substituting Eqs. (B4), (B5) into Eq. (A5) and taking 
into account that for a homogeneous system 02(77,77') = 
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G2(w — W r ), we arrive at the following result 



w = — 



2V9 



U\\ dU\ 



G 2 (€)de- (B6) 



Let us expand the pair correlation function, G 2 (£), m a 
Fourier series, 



G 2 (o = E e ^ G 2(k), 



(B7) 



and express W^ u , Eq. (B6), in terms of G 2 (k). 

First we note, that the following simple relation holds. 
Let F(|k|), where |k| = y/k^k^, be the Fourier compo- 
nent of a function i.e. 



F(\k\) = J e- <fc ^"F(|€|)d€. 



(B8) 



Then, the Fourier component of the function F(||^||) can 
be expressed in terms of F as follows 

/ e- ifc ^F(||€||)d€=-^F(||k||), (B9) 

where = y/g^^" (see Eq. (B5)), and 

||k|| = y/g^K (BIO) 

Substituting the expansion of Eq. (B7) into Eq. (B6) and 
using Eqs. (B8), (B9), we get the required representation 
for the interaction stress tensor 



•fl^dlkll 



G 2 (k). 



(Bll) 

In Eq. (Bll) the function w(|k|) is the Fourier compo- 
nent of the interaction potential, w(|£|), and w'{x) = 
dw{x)/dx. 

The stress tensor P^ u of Eq. (58) is the sum of 
Eq. (B3), and W^, Eq. (Bll). 



APPENDIX C: ELASTIC STRESS TENSOR IN 
COULOMB SYSTEMS 

The general expression for xc stress tensor, Eq. (110), 
can be represented in a much more simple form if the 
particles interact via Coulomb potential, w p = Ad/p d ~ 1 
(where A 3 = Ane 2 and A 2 — 2ire 2 ). Below we show that 
in this case the second term in Eq. (110) can be related 
to the potential energy, E pot , of a homogeneous electron 
gas with the density n/y/g. 

Let us represent the momentum integral in Eq. (110) 
as a sum of two terms 



where 



w„ v = w {1) + w (2) 



^=^f E^ q (v^^ 



(Cl) 



(C2) 



w$ = ~(d - i)f E A ^ G ? ( y/s^P^P fi }C3) 



To shorten the notations we retain only important mo- 
mentum dependence in the argument of the pair correla- 
tion function G^. 

Transformation of W$) , Eq. (C2), is straightforward. 
By changing the integration variables this equation can 
be reduced to the form 



= S 



9 

^ 2 



E 



A, 



-G?{p) (C4) 



Separation the integration over the modulus and the di- 
rection of momentum in Eq. (C4) yields the following 

result for Wfi) 



\g a pl a lPY 



E 



pot j 



(C5) 



where 1 is a unit vector (l 2 = 1), and the angle brackets, 
((...))i, stand for the averaging over the directions of 1. 

The momentum integral for wffi, Eq. (C3), can be 
reduced to a similar form. Let us first represent the de- 
formation tensor g^ in terms of its eigenvalues, A 2 , and 
eigen vectors, r]j^, 



_ \2 

g^v -^jVjuVjv 



(C6) 



Here j = 1, . . . , d labels the eigen vectors rjj that satisfy 
the completeness and the orthonormality conditions 



(C7) 



Tensor g^ has the same eigenvectors, while its eigen- 
values equal to 1/A 2 . Substituting the eigen vector ex- 
pansion of g^ v into Eq. (C3), and performing an obvious 
change of the integration variables we arrive at the fol- 

(2) 

lowing result for the tensor W^J 



(d-l)Af(ifel) 



[g a pl a lP] 



1 J±t~ ) ^p ot ' 



Combining Eq. (Cl), (C5) and (C8) we obtain the inter- 
action stress tensor W M „ in the following form 



W^ v = L llv (g a p)E potl 



(C9) 



where the calculation of the function L tlv (g a /}) involves 
only the angle integration. The angle integrals (factors 
with angle brackets in Eqs. (C5) and (C8)) are the scalar 
functions which depend only on eigen values of g^ v . For 
d = 2, 3 these integrals are reducible to a combination of 
the standard elliptic integrals. 
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